SINGULAR MODES OF THE ELECTROMAGNETIC FIELD* 
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Abstract. We show that the mode corresponding to the point of essential spectrum of the 
electromagnetic scattering operator is a vector-valued distribution representing the square root of 
the three-dimensional Dirac's delta function. An explicit expression for this singular mode in terms 
of the Weyl sequence is provided and analyzed. An essential resonance thus leads to a perfect 
localization (confinement) of the electromagnetic field, which in practice, however, may result in 
complete absorption. 
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1. Introduction. The ability to manipulate the spatial distribution of the elec- 
tromagnetic field is required in many practical applications. For example, one may 
wish to create an antenna with a very broad or a very narrow radiation pattern, 
localize and amplify light or transmit it along a predefined optical path, accelerate 
charged particles or keep them tight within a fusion chamber. In the presence of 
matter all possible spatial distributions are encoded in the spatial spectrum of the 
electromagnetic scattering operator. The most basic field distributions, which are 
easy to observe in microwave resonators, are called eigenmodes. Manipulation thus 
amounts to a clever excitation of a particular eigenmode or a combination of those. If 
almost all electromagnetic energy is carried by one of the modes, then we talk about 
a resonance. The concept of electromagnetic resonances and eigenmodes is a natural 
description of microwave resonators and waveguides 01 , microstrip lines |2| , and other 
simple, often, infinite, homogeneous, or periodic structures |S|, 01) El- In a recent 
paper [Sj we have generalized this idea for arbitrary dielectric objects of finite extent. 
The strongest point of our generalization was the incorporation of the full information 
about the spatial spectrum of the electromagnetic scattering operator [J] , which has 
both discrete eigenvalues and an essential (continuous) part. 

In some applications, e.g. optoelectronics, plasmonics, photonics, near-field op- 
tics, the electromagnetic field has to be confined within a very small volume of space, 
see e.g. HI, JO] . An observation made in |S] suggests that the modes associ- 
ated with the essential spectrum of the scattering operator may be highly localized 
in space. In particular, we argued that the resonance, which corresponds to the es- 
sential spectrum, coincides with the so-called plasmon observed at a plasma-dielectric 
interface. In research on metamaterials, where one strives for a negative permittivity 
material with vanishing losses, an unusual "anomalous localized resonance" is also 
encountered JT], ^21- Here we give a rigorous proof of the link between the local- 
ization or confinement of the electromagnetic field and the essential spectrum of the 
electromagnetic scattering operator. 



*This research was financed by the Netherlands Organization for Scientific Research (NWO) and 
by the Russian Foundation for Basic Research (RFBR). 

t Laboratory of Electromagnetic Research, Faculty of Electrical Engineering, Mathematics and 
Computer Science, Delft University of Technology, Mekelweg 4, 2628 CD, Delft, The Netherlands, 
n . V . budkoOtudelf t . nl 

^Department of Applied Mathematics, Moscow Institute of Radio Engineering, Electronics, and 
Automatics (MIREA), Verndasky av. 78, 117454, Moscow, Russia, absamokhinayandex.ru 



1 



2 



N. V. BUDKO AND A. B. SAMOKHIN 



Our previous work [7| was based on the Mikhlin's theory of singular integral 
operators ^^li and does not concern with the shape of the modes. In fact, there are 
very few examples of the analysis of the modes corresponding to the essential spectrum 
of operators. In quantum mechanics, where the spectral theory is very advanced, the 
discrete spectrum seems to be of more importance. There the essential spectrum 
is associated with the unbounded motion of particles and is, therefore, of little 
interest to physicists. In electromagnetics, however, it is the practical question of light 
confinement which is related to the essential spectrum, and the precise mathematical 
description of the associated modes is rather important. 

To recover the structure of the modes corresponding to the electromagnetic es- 
sential spectrum we resort here to the Weyl's definition of spectrum, see e.g. [T^ . 
|15|. This definition states that a number A is in the spectrum of operator A if and 
only if there exists a sequence {vl/n} in the space X such that 

ll*n||=l, (1-1) 

and 

lim ||A^'„- A*„|| =0. (1.2) 

n — ^CJO 

Furthermore, A is in the essential spectrum, if there is a singular sequence satisfying 
()l.l|l - (|1.2|l . i.e. a sequence which contains no convergent subsequence. If X is a 
complete Hilbert space, then, obviously, such singular sequence will not (strongly) 
converge to any function from X, although, it may weakly converge to zero. As 
we know, sequences that do (strongly) converge to some function on X generate 
eigenf unctions or eigenmodes corresponding to the point spectrum - eigenvalues. By 
analogy we may associate essential modes with the essential spectrum. An alternative 
term, which, perhaps, better reflects the nature and structure of the particular modes 
obtained here, would be singular modes. 

2. The square root of the delta function. Before going into the details of 
the electromagnetic case we shall introduce a peculiar function that is later used to 
generate the essential modes. On one hand, the sequence of such functions should be 
singular, as required by the definition mentioned above. On the other hand, appli- 
cation of (|1.2|l to electromagnetics (in L2-norm) and consistency with the previously 
obtained results require that the square of this function should behave like the 
Dirac delta function, i.e., must have the sifting property. Hence, what we need is a 
square root of the delta function. 

From time to time the square roots of delta functions appear in literature. Mostly, 
though, just as a curious example of a non-convergent sequence, see e.g. jltij (p. 299) 
and .l7| (p. 81). A more modern and rigorous approach to such functions is the 
Colombeau algebra m| , where one studies the products of distributions and encoun- 
ters m-singular delta functions, which are almost identical to what we are after. At 
present, the main applications of the Colombeau algebra are: nonlinear equations 
(e.g. hydrodynamics, elastodynamics and general relativity), singular shock waves in 
nonlinear conservation laws, and propagation of delta-like waves in linear media with 
discontinuous parameters. In other words, this algebra is applied whenever a poten- 
tially meaningless product of generalized functions is stumbled upon. The present 
apparent emergence of the Colombeau algebra in the recovery of the essential spec- 
trum seems to be new and can, probably, be generalized. This paper, however, does 
not attempt such a generalization, nor should it be considered an introduction to the 
Colombeau algebra. 
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Unfortunately, we could not find any explicit derivation of the required distribu- 
tion in the literature. The available one- dimensional and scalar three-dimensional [Tl) 
(pp. 74-75) examples are of no use to us, since the electromagnetic essential spectrum 
is a purely three-dimensional phenomenon and our function and its Fourier transform 
had to have a very special vectorial structure. In view of the potential usefulness of 
the obtained result in other areas of research, we have decided to devote this entire 
section to the analysis of the square root of the three-dimensional delta function. The 
proof of the following theorem is, however, rather technical and can be skipped in the 
first reading. 

Theorem 2.1. The vector-valued function 



^'(a,x,Xc 



1/2 



7r-3/4a5/4 



(x - Xc ) exp ^ |x - Xc 



(2.1) 



where x, Xc G M.^ and a > 0, has the following properties: 

1. It is normalized in the sense that 

||*(a,X,Xe)|l2 = 1- 

2. The sequence of such functions 

5'(q!„,x,Xc), q;„ > q;„_i, n=l,2, 

does not have a convergent subsequence. 

3. Its Fourier transform is given by 



(2.2) 



(2.3) 



*(a,k,Xc 



'2 



1/2 



^_3/4^-5/4i^g^p(_^|j^ 



1 



2a' 



ik • Xr 



(2.4) 



4-. It is a generator of the square-root of the Dirac delta-function, i.e., with any 
bounded continuous function /(x) its square has the sifting property 

lim / /(x)|M/(a,x,Xe)|' dx = /(xe). (2.5) 

5. It is orthogonal to bounded vector-valued functions, i.e., for any |V(x)| < oo, 

X e R^, 

lim(V,*)=lim / *'^(a,x,Xc)V(x)dx = 0. (2.6) 

xGR3 

6. It is 'invisible' to weakly singular operators with finite spatial support, i.e. 



lim 

a — *oo 



]K(x,x') 



c'l/3 



^'(q;, X, Xc) dx 



(2.7) 



where K(x, x'), x, x' G is a bounded tensor-valued function, /3 < 3, and 
the norm is defined over the spatial support D. 
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Proof. 

1. To prove the normalization property we simply compute 

||*(a,x,Xc)||^ = ^7r-3/V/2 J |x-xepe-«l---=l'dx 

xeiR3 

oo 7r 2-77 

= |7r-3/2a5/2 J J J r^e-"'-' sinM^dMr (2-8) 
where we have used the following standard integral: 



r=0 6»=0 V=0 

2 -3/2 5/2. (4-1)!! /^y/^ . 



OO 



2. Suppose that there is a subsequence \l/(a„, x, Xc), where a„ > a-n-i, n = 
1,2,..., which converges in norm. Then, for any e > there exists A'' such that for 
all m,n > N we have 



|*(a„,x,Xc) - *(a„,x,Xc)f < e. (2.10) 



However, in our case we obtain 



|^'(a„,x,Xc) - *(a„,x,Xc)||^ = J |*(a„, x, x^) x, x^)]^ dx 



= j |*(Q!m,x,Xc)|^ dx+ j |*(a„,x,Xc)|^ dx 

xeR3 xeM3 

-2 j *^(aTO,x,Xc)*(a„,x,Xc)dx 



(2.11) 

eM3 



5/4 5/4 
2 2 CKn 



xeK^ 



y |*(am + Q!„,X,Xc)|^ dx 



5/4 5/4 

2-2- "™ 



{am + OLnfl"^ ' 



If we now fix n > iV and let m > iV go to infinity, then the last term tends to zero 
leaving us with a constant (two), which, obviously, cannot be made smaller than an 
arbitrary e > 0. This proves that there are no convergent subsequences. 
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^(ajk.Xc) = (27r) j *(a, x, Xc) exp(-ik • x) dx 

x£M3 

= (27r)-3/2 (J^ TT-^/^a^/* y (x - Xe) exp (-||x - x^l^ - ik • x) dx 

xeiR3 

7r-3/V/*expHk.Xe) J yexp(-||y|2-ik.y) dy 

yeK3 

7r-3/V/4expHk.Xc)(zVk) j exp (-||y|2 - ik ■ y) dy 

oo 

7r~^/'^a^/'^exp(-ik-Xc)(iVk) j exp (^-^yf - ikiyij dyi 



= {2n) 
= {2n) 



-3/2 / f 



-3/2 / f 



1/2 



1/2 



1/2 



cx> oo 

/ (~f ^2 - ^^2^2) dy2 X exp (-^2/3 - ^^sj/s) dya. 



(2.12) 



Each of the one-dimensional integrals above gives 





' f 2 


.2 Y 


exp 











dy„ 



00 

/ 



exp 



dyn 



= exp 



= exp 



00 

-i*") / 

— CXD 



exp 



exp 



~ TT I 2/n H fen 

2 \ a 



dyn 



277 ' V2 



ira 



oo+ib 

?) ^''^("i^") / exp [-7r(^ + d(^ + i6) 



-oo+i6 



1 



(2.13) 
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Substituting (|2.13|l in (|2.12|) we arrive at our result 
*(a,k,Xc) = 

i(2^)-3/2 Q\ ^-3/4^5/4 exp(_ik • Xc)Vk 



-i(27r) 



0\ 1/2 /97r\3/2 1 

-3/2ff'l ^-3/4^5/4 / \ i 



1/2 



exp I - ik • Xc 



Note that with this choice of the Fourier transform we also have 

*(a,k,Xc) 



(2.14) 



(2.15) 



4. To prove the sifting property we split the integration domain in two parts, i.e., 
lim / /(x)|*(a,x,Xc)|' dx 

Q — ►OO J 



lim 

a. — >OQ 



J /(x)|vl/(a,x,xe)|' dx+ J /(x)|vl/(a,x,xe)|' dx 

xeR3\y(5) xev{S) 

lim / /(x)|*(a,x,Xe)|' dx+ lim / /(x) |*(a, x, Xe)|' dx, 



<:GR3\y(5) 



ieV{5) 



(2.16) 



where V{S) is some volume surrounding the point Xc. Now we shall choose V{d) in 
such a way that the first (outer) term in the expression above gives zero. We start by 
considering a general case where x S V{6), if |x — xd < 6, and S is some function of 
a. Then, 



lim / /(x) |*(a,x,Xc)|^ dx 

a — *oo J 

x6R3\V(5) 

<max|/(x)| lim / |4'(q!, x, Xc)|^ dx 

xGR3 a^oo J 



<;eR3\y((5) 



max I /(x) I lim -n-^/^a^^'^ 

xgR3 a^oo 3 



(2.17) 



r=s e=a 1^=0 

OO 



— 7r-3/2niax|/(x)| lim a'/^ f 

r—6 



dr. 



Successive integration by parts gives 



r'e "" dr = — ^ / e 
4a2 

r—5 r—S 



^ , <53 3(5 
'^''+'2^ + 4^'^ 



(2.18) 
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Using this result we continue to analyze the upper bound of (|2.17|l as follows: 



lim 

a — >oo 



|*(a,X,Xe)|' dx= lim ^^-3/2^5/2 / ^4g-ar-j^ 
Of — >oo 3 



tGM3\y((5) 



r—S 



^^-3/2 lini a5/2 

3 a — >^oo 



3 

4q:2 



e""'' dr + — 



6^ 35 



2a 4a2 



-aS-' 



r—S 



(2.19) 



= ^ lim a^/^ / e^"''' dr + lim ( 

r—S 

Now, choosing, for example, 5 = a^^/^, we arrive at 
lim / X, Xc)|^ dx 

a— ♦oo J 



2 4 



xeR3\y(5) 
2 



L lim / e-^^^^" d(rVa) + lim f + ^a'A e""' 

TT a^oo J a^co 3y7r \2 4 / 



oo 

2 . /■ _ 
< — ^ lim / e ^ dz 



lim e-" = 0. 



(2.20) 



Hence, with this particular choice of V{S) the first (outer) term in H2.16(l is zero. Now 
we shall use the same V{5) in the second (inner) term. Taking into account that /(x) 
is a continuous function, and that with our choice of d the integration volume V{S) 
tends to the point Xc, we can apply the mean-value theorem, i.e.. 



lim / /(x)|*(a,x,Xc)|' dx= lim /(x„) / |*(a, x, Xc)|' dx 

a — ♦oo / a — ♦oo / 



iev{S) 



xev{S) 



/(xc) lim / |*(a,x,Xe)|' dx, 

a^oo J 



(2.21) 



xev{S) 



where Xq, e V{6), and x^ — )■ Xc as a ^ oo. Thus, to prove the sifting property it 
remains to show that 



-1/3 



lim 



/ |*(a,x,xc)|^ dx= lim a^/^ / r^'e""'^' dr 
J 3V7r "^oo J 

-V(S) r=0 



xev{S) 



,1/6 



,1/6 



lim / (rVa)^e-(''^)' d(rV^) = ^ lim / z^e'^" dz (2-22) 



3y/Tr a- 



z=0 
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Finally, we remark that the above proof holds with any 5 = a~^/™, where m is an 
integer m > 3. 

5. We prove the orthogonality property by considering the absolute value of the 
dot-product, i.e., 

lim (V(x),*(a,x,Xc)) 



OO TT 27r 



lim Jt^-^/^a'/^ 



r=Oe=0 (p=0 



@'^{$,ip)V{r, 6,1^^6-^'' sinMiysdMr 

OO 

< Cmax|0^(6',v5)V(x)| lim a^^^ [ r^e'^'^'dr 



(2.23) 



r=0 



= Cmax|0T(0,^)V(x)| lim a^/*— 1— 
= 2Cmax|e'^(6»,</j)V(x)| lim a-3/4 ^ 

xGM^ a— >(X) 

where < C < oo, = (x — Xc)/|x — Xc|, and the following standard integral was 
used: 



oc 

/ 



^2„+ig-pr^d^^__^ P>0, n = 0,l,2,. 



2p 



(2.24) 



6. In the case of a weakly singular integral operator with finite spatial support D 
we proceed as follows: 

2 



lim 

a— ^cxj 



f K(x,x') ^ , , _^ 

J Ix-x^ ^"''''^^^^'' 



< 



xe£> 



x'eD xeu 



x'eD xeD 



x'£r> xgd 



]K(x 


,x') 


|x- 


x'f 


K(x 


,x') 


|x- 


x'|/5 


K(x 


,x') 


|x- 


x'|/3 


K(x 


x') 


|x- 


x'|/5 



dx' 



dx' (2.25) 



dx' 



dx', 

J Ix-x'l^^ ■ ■ " 

x'eD x£D 

where the original function ^ is split into two complementary parts with respect to a 
small volume around Xc in such a way that: 

'i'(a,x,Xc) = ^'i(q:,x,Xc) + *2(a,x,Xc), 

*i(a,x,xc) =0, X e R3 \^ t/((5), (2.26) 
*2(q;,x,Xc) = 0, xeV{S). 
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The last integral in (|2.25|l is estimated like this 

f f K(x^ 



x'eD 



dx' 



(2.27) 



lim \\K^2f < \\Kf lim ||*2||' = 0, 

a — ^oo a — ^oo 



where we have used the fact that the norm of a weakly singular operator on D is 
bounded, and the previously derived property 12.19(1 - 12.2011 . This means that we 
take > 3. The remaining integral in H2.25|l requires considerably more 

work. We shall split the domain of integration over x' in two parts using yet another 
small volume V{S') surrounding the point X;.. Then, taking H2.25|l into account, we 
obtain 



f f IK x,x' ^ , 
hm / / --T*i(a,x,Xc)dx 

a^oo J J |x-x'|P 



dx' 



lim 



x'eD\V{5') 



K(x,x') , 

— — -^^'i(a,x,Xc)dx 



dx' 



(2.28) 



lim 



IK(x,x') . . ^ 

— -— ^*i(a,x,Xe)dx 



cey(<5) 



dx' 



Proceeding with the first of the above integrals we apply the Caushy-Schwartz in- 
equality and arrive at 



lim 

a — *oo 

x'eD\V(5') 

< lim 



IK(x,x') 



-g*i(a,x,Xc)dx 



xey(5) 



dx' 



|K(x,x')0|' 



j--^dx I |^'i(a,x,Xc)|^ dxdx' 

a^oo I I IX — X ' 



x'eD\v{S') xev(5) 
< max |K(x,x')0|^ 



X lim 



x'eD\v{S') xev(S) 



iev{5) 



(2.29) 



Now, if we choose the spherical volume V{6') with radius 6' to be larger than the 
spherical volume V{S) with radius S, then for x S V{6) and k' G D \ V{d'), the 
distance factor |x — x'| will be bounded from below by the difference of the radii of 
the two volumes. In the estimate (|2.27f) we have used S — a"^/™, m > 3. Hence, for 
sufficiently large a we can choose S' — (5^/" ~ a^^/inm) ^^^^ integer n > 2. In this 
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case we have 



lim 

a — >oo 



< lim 



/ / |x-x^|2/3 / l*i("'X'Xc)|^ dxdx' 

(<5) 

y" J dxdx' y" |*i(a,x,xc)|' dx 



x'eD\y(<5') xev(S) 



xev{S) 



max 



a^coxeV{S),x'eD\V{S') |x — x' 



< C lim 



C lim 



C lim 



<;'er>\y(<5')xe\/(<5) 

^3 



^(l-3n)/(rim) 



2/3 



= C lim 



xey(<5) 

(53-l/n 



(l-51-l/n) 



2/3 



a^oo _ Q,(l-n)/(nm)~) 



2/3 



= 0, 



(2.30) 



i.e., the first of the two integrals in (|2.28|) is zero. Applying the Caushy-Schwartz 
inequality we estimate the last integral in H2.28|l as follows: 



lim 



x'ei^(5') 



— -— ^*i(a,x,Xe)dx 



xev(S) 



dx' 



x'ev{6') xeviS) 

< L lim {6')^-^ 

a — >oo 

x'ey(<5') xev(i5) 

< M lim (J')^^^-'^) / |*i(a,x,Xe)|^dx 



|^'i(a,x,xc)|2 



dxdx' 



(2.31) 



xev(<5) 

= N lim (y)6-2/3 = N lim a-(6-2/3)/(nm) ^ 0^ 

a — >oo a — >oo 

where 6' = a~"'^/*^"'"\ m > 3, n > 2, while < /3 < 3 by the conditions of the theorem. 
Thus, we have shown that H2.25|l is, indeed, zero. □ 

3. Electromagnetic singular modes. Consider the volume integral equation 
of electromagnetic scattering on a nonmagnetic object of finite spatial extent D: 

E'''(x,w) = E(x,w)- 

[ft2(u;) + VV-] J g(x-x',c.)x(x',c.)E(x',c.)dF, (^.1) 
x'eD 

where E'" and E are the incident and total electric fields, correspondingly. This 
equation is obtained directly from the frequency-domain Maxwell's equations and 
takes into account the radiation condition at infinity in a most natural form. The 
medium parameters are contained in the contrast function x, which in terms of the 
complex permittivity function e will look like 



X X, a; = -^-^ - 1 = er{x,uj) - 1. 
£0 



(3.2) 
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The vacuum wavenumber is ko — lo/c, and the scalar Green's function is given by 

gifco|x| 



?(x,a;) 



47r X 



(3.3) 



Carrying out the two spatial derivatives we arrive at the following singular integral 
equation: 



E-(x,c.) 



1 , 

1 + gXlX,^ 



E(x,w) 

lim / (Go(x-x')x(x',w)E(x',w) dx' 

(S-.0 J 

x'en\|x-x'|<<5 
J Gi(x-x',t^)x(x',w)E(x',tj) dx'. 



(3.4) 



x'eD 

For the definitions of the Green tensors Go and Gi we refer to our previous publication 
on this subject |7|. 

Theorem 3.1. The vector-valued function ^'(a,x, Xc) defined in Theorem \2.1\ 
generates the essential mode of the electromagnetic field corresponding to the point of 
essential spectrum 



Acss — £r(Xc, uj). 



(3.5) 



Proof. Since we have already established the normalization (Property 1) and the 
singularity of the sequence (Property 2) in Theorem 12.11 we only need to prove the 
following analogue of (|1.2|) : 



V= lim 

a — !-oo 



1 + 3X(x,w) 



*(a,x,Xc) 

lim / Go(x — x')x(x', (jj)'I'(q, x', Xc) dx' 

5^0 J 



x'eD\|x-x'|<(5 



j Gi(x- x', w)x(x',cj)*(q!,x',Xc) dx' - A*(q!,x, Xc) 



x'eD 
= lim 



D 



3 + 3S.(x)-A 



^'(a,x,Xc) 

-lim / Go(x- x') [er(x') - er(x) +er(x) - 1] «'(a,x',Xc) dx' 

<5^0 J 



x'eD\v{S) 



J Gi(x- x',cj) [£i.(x') - 1] *(a,x',Xc) dx' 



x'eD 



= 0, 



D 



(3.6) 



12 



N. V. BUDKO AND A. B. SAMOKHIN 



for A — ei-(xc,w). The Li norm is taken over the finite spatial support D. First, we 
rearrange and decompose it into separate terms 

V < lim ||[e,(x)-A]vI/(a,x,Xc)||^ 



lim 



+ lim 

a— ^oo 



[1 - 



^4'(a,x,Xc) + lim / (Go(x - x')1'(a, x', Xc) dx' 

6 <5^0 J 



'(iD\V{S) 



D 



lim / Go(x- x') [e(x,w) - e(x',w)] «'(a,x',Xc) dx' 

x'eD\v{S} 

J Gi(x - x', uj) [e,(x') - 1] ^{a, x', Xe) dx' 



x'eD 



D 



(3.7) 



With £i.(x, Holder-continuous in M"^ all integral operators in the last term are 
weakly singular. Hence, from the Property 6 of Theorem 12. II the last term in (|3.7|l is 
zero. 

From the second term in (|3.7|1 we only consider the expression in the large square 
brackets, which we extend to R'^. Denoting by JF{. . . } and !F~^ {■ ■ ■} the forward 
and inverse three-dimensional Fourier transforms, we obtain 

^^'(a,x,Xc) + lim / (Go(x - x')^'(a, x', Xc) dx' 

6 S^Q J 

x'GR3\y(5) 

^J'^{*(a,x,Xe)}+J^i Inn J Go(x - x')^'(a, x', Xc) dx' 



|^'(a,k,Xc) 



.F-M#(a,k,Xe) 



kkT ^ 



x'eR3\y(5) 
*(a,k,Xc) 



*(a,k,Xe) ^ =0, 



(3.8) 



where we have used the previously derived result about the Fourier transform of a 
simple singular operator [7], the explicit form of 4'(a, k, Xc) (Property 3), and the fact 
that k"'"k = |kp. Hence, the second term in H3.7|l is also zero. 

Finally, applying Property 4 (sifting) of Theorem 12. II we see that the first term 
in H3.7|l is zero, if equality H3.5|l holds. □ 

4. Conclusions. Now we know that the electromagnetic essential (singular) 
mode is, in fact, the square root of the delta function. But what does it mean from 
the physical point of view? Is it possible to have an "essential" resonance? Under 
which conditions? Can the electromagnetic field be confined to what seems to be a 
single point in space? These are open questions. However, one important conclusion 
can be deduced from the basic principles of the electromagnetic theory. 

The resonant excitation of a singular mode is only possible if the essential spec- 
trum contains zero. Physically this means that er(xc,a;) = at some point Xc in 
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the scattering domain D. It may seem improbable, but one has to remember that 
the dielectric permittivity is an effective macroscopic parameter, which has very little 
meaning for the microscopic induced current densities at the atomic level. Moreover, 
in the classical Lorentz theory of atoms, the real part of the dielectric permittivity 
in a dispersive medium can vanish and even become negative (the so-called Lorentz 
absorption line). Although, in that case the Kramers-Kronig relations predict an in- 
crease in the imaginary part of macroscopic er, which corresponds to the absorption 
of the electromagnetic field. Turning it all around we can say that even if at a micro- 
scopic level we could have ei-(xc,w) = and excite the corresponding singular mode, 
it should be absorbed to preserve the macroscopic Kramers-Kronig relations. 

In principle, the very fact that singular modes are outside the Hilbert space, where 
all 'proper' solutions of the Maxwell equations live, tells us that the excitation of sin- 
gular modes can be considered as some kind of 'deflation' of the electromagnetic field. 
Recall that, due to the normalization Property 1 of Theorem l2.1l the electromagnetic 
energy associated with a singular mode is well defined. 

It is also interesting to note the direct relation of the electromagnetic essential 
spectrum and its singular modes to the pseudospectrum and wave-packet pseudo- 
modes [21 • There is an obvious similarity of the Weyl's definition (|1.2I) and the 
definition of the pseudospectrum, where instead of zero one should simply put a small 
e in the right-hand side of (|1.2|1 . Subsequently, we arrive at two distinct possibili- 
ties. The first is where n — > (X), i.e., in our case a oo. Then, points Aps satisfying 
I Aps — er(x, i-u)\ < e will belong to the pseudospectrum, while the corresponding modes 
will be singular. The second case is where Aps — Xcss — £r(x, a;) or very close to it, 
but a < 5{e). In this case, we stop the sequence of 4'(a,Xc,x), at some finite a, 
for which the norm in H3.6(l equals e. Although, it is difficult to derive an explicit 
relation for (5(e), we can anticipate that ^'((5(e), Xc, x) will be highly localized in space 
around the point Xc. In this case the mode is not singular and belongs to the Hilbert 
space. These two physically distinct possibilities emphasize the nonunique nature of 
the pseudospectrum as it is defined in |19| and elsewhere. 
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